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Abstract 

We present a molecular dynamics study of reentrant nematic phases using the Gay-Berne-Kihara 
model of a liquid crystal in nanoconfinement. At densities above those characteristic of smectic A 
phases, reentrant nematic phases form that are characterized by a large value of the nematic order 
parameter S ~ 1. Along the nematic director these "super nematic" phases exhibit a remarkably 
high self-diffusivity which exceeds that for ordinary, lower-density nematic phases by an order of 
magnitude. Enhancement of self-diffusivity is attributed to a decrease of rotational configurational 
entropy in confinement. Recent developments in the pulsed field gradient NMR technique are 
shown to provide favorable conditions for an experimental confirmation of our simulations. 



PACS numbers: 61.30.Hn,66.10.C-,64.70.M- 61.30.Gd 



In the context of phase transitions "reentrancy" refers to a nonmonotonic variation of 
an order parameter with the thermodynamic field driving the transition. Reentrancy is 
ubiquitous in the physics of thermal many-particle systems. It arises under quite disparate 



physical ^conditions encountered, for example, in quantum gases two-dimensional charged 



colloids {2], or relativistic scalar field models [3j. As far as soft matter is concerned reentrancy 



has been reported for self-assembled supramolecular structures 



oleophilic surfaces jfj], and novel discotic and calamitic liquid crystals [7j|. In fact, since the 
first observation of reentrant nematic (RN) phases in a seminal paper by Cladis 8[ reentrancy 
in liquid crystals seems to have received most of the attention. This is most likely because 
of the abundance of phases exhibited by these materials. For example, reentrant phase 
;ransitions have been reported for the isotropic phase of mixtures of discotic liquid crystals 
9], the ferroelectric transition in syn- and anticlinic smectic C phases [loj], the cholesteric- 
to-blue phase transition in chiral liquid crystals 11], and for nematic (N) phases 12]. 

Despite the variety of systems and thermodynamic conditions under which reentrancy 
in liquid- crystalline materials arises comparatively little attention has been paid to the 
dynamics of reentrant phase transitions. For example, distinct differences in the molecular 
dynamics in the N and RN phases can be concluded from corresponding changes in the 
nuclear magnetic relaxation times reported in Refs. I31 15. 

Whereas most earlier work on reentrancy of phase transitions in liquid crystals is experi- 
mental in nature comparatively little attention has been paid to this fascinating phenomenon 
from a theorist's point of view. The most recent theoretical study employs isothermal- 
isobaric and canonical ensemble Monte Carlo (MC) simulations to investigate the nature of 
the smectic A(smA)-RN phase transition for a bulk system of hard ellipsoids with square- well 
attraction 16[ where earlier theoretical studies are briefly reviewed, too. Unfortunately, the 
model employed in Ref.[l6|is somewhat artificial in assuming that the ellipsodal molecules are 
always oriented in a perfectly parallel fashion such that all rotational degrees of freedom are 
always "frozen" irrespective of the thermodynamic conditions. Therefore, this study seems 
only of limited use to elucidate properties of RN phases at a molecular level. Moreover, the 
authors do not consider dynamic features of RN phases. 

Therefore, we show here that the smA-RN phase transition causes a dramatic increase 
in the self- diffusion of the molecules in the direction of the nematic director n. Our model 
system consists of soft spherocylinders interacting via the so-called Gay-Berne-Kihara pair 



wetting phenomena at 
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2 



potential [17]; specifically, we use the GBK(6,5,2,1) version of the model in the notation of 
Ref. [l?! Data will be presented for a system of spherocylinders confined to a slit-pore with 
structureless walls, separated by a distance s z along the z-axis. The fluid-solid interaction 



is described by the surface- averaged potential 



Uf s = 4e fs p s cr 



a 



10 



a 



(i) 



In Eq. (CQ), ef s determines the depth of the attractive well, p s cr 3 = 2- 1 / 3 is the areal density 
of a layer of substrate atoms, d w is the minimum distance of a spherocylinder from either 
substrate, and g (u) is the so-called "anchoring function" where the unit vector u specifies the 
orientation of a molecule. The anchoring function is introduced to discriminate energetically 
specific orientations of molecules with respect to the substrate plane following a suggestion 
of Steuer et al. 



191 ] . Here, g (u 



u; 



such that an orientation parallel with the 



substrate plane is favored 20]; throughout this work s z = 19a is used. 

We employ both extensive isothermal-isobaric MC and microcanonical molecular dynam- 
ics (MD) simulations to locate the smA-RN phase transition and to study changes in mass 
transport accompanying that transition. The main purpose of employing MC is to provide 
suitably equilibrated starting configurations for subsequent MD simulations. In addition, 
MC is used to independently verify the correctness of the MD simulations through a com- 
parison of equilibrium properties obtained in both types of simulations. In MC we use a 



standard algorithm 21] but allow the side lengths of the computational cell to vary inde- 



pendently to preserve the in-plane isotropy of the pressure tensor even in highly ordered 
rases; MD simulations are based on an implementation of the velocity Verlet 



confined p 
algorithm 

perature T, and (<r 2 m/e) 1//2 for 



221 ] . We employ the customary dimensionless units of a for length, e/k-Q for tem- 



;ime t, where m — 1 is the spherocylinder mass. Values for 
a and e are taken from Ref. [17|. Our simulations comprise N = 1000 molecules of length 
L = 6. Interactions are truncated beyond a minimum distance d m = 3. In MD the starting 
configuration is further equilibrated using a velocity rescaling scheme to the desired T at 
constant volume. Throughout this work we employ a time step of At = 10 -4 to integrate 
the equations of motion. 

To characterize the degree of nematic order we follow established protocols and compute 
the nematic order parameter S as the ensemble average of the largest eigenvalue of the 



3 



O.Od 



o.nx 



0.1 0.12 0.14 

P 



FIG. 1: (Color online) Plots of S (O,*) and A (□,■) as functions of number density p for confined 
liquid crystals. Open and filled symbols refer to T = 4.0 and 6.0, respectively. 




FIG. 2: (Color online) Configuration "snapshots" from MD; p = 0.110 (left, N), 0.133 (middle, 
smA), 0.139 (right, RN), T = 4.0. Color code: Ui - n = , «j ■ n = 1 ,i = l,...,N. 



instantaneous alignment tensor 23 j 

Q = 




(2) 



where "(g)" denotes the direct product and 1 is the unit tensor. Hence, Q is real, symmetric, 
and traceless. Its eigenvector associated with the largest eigenvalue corresponds to n. The 
layering characteristic of smA phases is quantitatively described through the function [lj| 



X(d) 



1 

N 



N 



J^exp 



3=1 



2ni [Tj ■ n) 

d 
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where d is the spacing between layers and the smectic order parameter A is defined as the 
maximum of A (d) in the interval [L — (,L + Q where ( = 0.05. Focusing on T = 4.0 first, 
plots in Fig. [TJ show that S is rather small up to a number density p < 0.09 indicating 
that the confined fluid is in its isotropic phase. Beyond p > 0.09, S rises steeply assuming 
values characteristic of nematic phases (S > 0.4) and levels off as p increases. At the highest 
densities considered S ~ 1 which reflects a nearly perfect orientation of molecules. Over 
the same range of densities, A is very small up to values of p ~ 0.115, indicating that the 




FIG. 4: D\\ [see Eq. ©] as a function of p for the confined liquid crystal; (O) T = 4.0, (•) T = 6.0. 
Inset is a magnification of the T = 4.0 isotherm around the smA-RN transition. 

confined fluid does not form any smectic layers. As the density increases further smectic 
layers are forming revealed by A > 0.5. Interestingly, as the density keeps increasing beyond 
p ~ 0.135, A drops to a low value of about 0.1 whereas S increases further towards its 
maximum value of about 1. Hence, p ~ 0.135 demarcates the formation of a RN phase. 
The structural characteristics of confined phases pertaining to the stable regime of N, smA, 
and RN phases is illustrated by "snapshots" of configurations generated in MD (see Fig. [2]). 
Notice the presence of distinct molecular layers in the smA state which disappear once the 
RN state forms. Comparing snapshots for typical N and RN phases it is apparent from 
Fig. [2] that the latter exhibit much more orientational but roughly the same positional order 
compared with the former. 

The key result of our study concerns the self- diffusion of molecules in the direction of 
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their orientation. To that end we introduce the parallel mean-square displacement (MSD) 

The angular brackets indicate an average over the iV molecules and a sufficiently large 
number of time origins referred to by subscript "t" . In Fig. [3] we present plots of the 
MSD for typical N, smA, and RN phases where different time regimes can be identified. 
In the double-logarithmic representation short-time slopes of the MSD's exceed their long- 
time counterparts on account of the initial ballistic motion of molecules in all three phases. 
As expected, the MSD's are indistinguishable for t < 0.1. In the limit of large r the time 
dependence of the MSD's corresponds to diffusive motion, that is the MSD's depend linearly 
on t. The MSD for the smA phase exhibits a plateau at intermediate times which eventually 
gives way to diffusive motion at long times. The plateau reflects the presence of smectic 
layers that hinder molecular motion in the direction perpendicular to the plane of adjacent 
smectic layers (see Fig. [2]). In the limit r — > oo reliable values of D» can be extracted from 



the plots in Fig. |3] via the Einstein relation 24 
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lim^(ArJ(r)> t (5) 

A comparison of plots in Figs. [1] and H] reveals that in the N phase, D\\ is small but nonzero. 
It attains a nearly vanishing value in the smectic phase on account of layer formation (se e 
Fig. [2]) which blocks mass transport efficiently in the direction of the layer normal 24]. 
Most importantly, however, compared with the smA phases D\\ increases dramatically by 
several orders of magnitude when the RN phases become thermodynamically stable (i.e., 
f° r P ^ 0.135, see inset in Fig. H]). Because of the unusually large D\\ in conjunction 
with high nematic order (S ~ 1, see Fig. [1]) we call these high-density nematic phases as 
"supernematic" . One also notices by comparing data for different T in Fig. H]that both data 
sets show a similar large increase of D\\ beyond a certain density threshold. However, the plot 
for T = 6.0 exhibits nonzero values of D\\ over a density range where the corresponding curve 
for the lower T = 4.0 drops to zero. This is due to the absence of smA phases at T = 6.0 (see 
Fig. [1]). Interestingly, the experimentally determined phase diagram presented by Guillon 
et al. also suggests absence of intermittent smA phases and a continuous transition from N 
to RN phases under suitable thermodynamic conditions 



. Nevertheless, the increase in 
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FIG. 5: D\\ as a function of o^ 1 for T = 4.0 (O) and T = 6.0 (•). The linear fit ( — ) is rationalized 
in the text. 

orientational order for p > 0.135 and T = 6.0 causes D\\ to increase equally strongly. Hence, 
at sufficiently high T one may observe supernematic features without reentrancy. 

The apparent dramatic increase of self-diffusivity in the RN phase can be attributed to a 
loss of rotational configurational entropy S TC due to both increase in density and presence of 
solid surfaces. We rationalize this by assuming a characteristic time interval t\\ associated 
with the onset of diffusive motion in the direction of n such that D\\ oc 1/tii. For r > tii 
the probability V that a molecule has traveled a distance Arii from its origin at r = in 
the direction of n should then also be inversely proportional to tii. Intuitively, one expects 
V (Ary) to be larger if the alignment of molecules with n is more pronounced on average, 
that is the larger S is. However, a larger value of S implies a lower rotational configurational 
entropy S IC such that V oc exp (— S rc /k B ) using standard statistical-physical reasoning \^\ . 
This then suggests D\\ oc exp (— S tc /k B ). We estimate S rc via 

S IC = -k B J &6 P{6) In P(9) (6) 

where P(9) is the distribution of angles cos 6^ = ui ■ n. Assuming P(6) to be Gaussian 
with a standard deviation of a rc it is easy to verify that <S rc oc In a r r . A similar relation was 



28l assuming a Gaussian 



obtained for the configurational entropy of a macromolecule in Ref. 
distribution of relevant internal degrees of freedom. Hence, the above line of arguments 
suggests D\\ oc a" 1 which is supported by plots in Fig. [5j 

In summary, we have shown for the first time that supernematic liquid-crystalline phases 
exhibit unusually large self-diffusivity which can be explained in terms of a substantially 
reduced S rc . Here the suppression of S TC is assisted by confinement to a nanoscopic slit- 



pore where the solid substrates favor a planar arrangement of the molecules with respect 
to the substrate plane. Because the solid surfaces of the slit-pore may be viewed as the 
representation of an external field superimposed onto the intermolecular interactions, we 
anticipate the results presented here to be generic in that they should persist in other fluids 
composed of anisometric molecules that are exposed to external fields (e.g., dipolar fluids in 
external magnetic fields). Hence, we believe our results to be important for a broad range 
of liquid- crystal applications ranging from lubricants in nanotribology 29| over nanosensors 
30( to photonic 3l| and organic electronic devices 32| where the mobility of molecules plays 
a key role. 

Though the experimental data so far communicated in the literature do not yet provide 
a rigorous confirmation of our theoretical predictions, some of these data are in reasonable 
qualitative agreement with our simulations. For example, extrapolating longitudinal relax- 
ation rates reported in Ref. [14j from the RN to the N phase yields values markedly below 
those in the N phase. If referred to the same T these relaxation rates correspond to corre- 
lation times that are notably shorter in the RN compared with the N phase. However, this 
general interpretation of NMR data remains speculative as long as translational diffusion 
has not definitely be identified as the process governing the observed relaxation. 

Direct evidence for the relation between translational diffusion and NMR data can be pro- 
vided by the pulsed field gradient NMR (PFG NMR) technique 33| which records molecular 
displacements typically over a /im range. For example, PFG NMR has been applied to di- 
rectly assess the diffusion tensor upon entering the N phase [34j. In these studies diffusion 
in the direction of the molecules' long axes was found to increase with increasing nematic 
order. By the same technique the diffusivity of n-alkanes in nanochannels was found to 



increase with increasing orientational order 



351 ] . These findings are in line with our data 



where enhanced molecular ordering is accompanied by increasing diffusivities in the direc- 
tion of n. The powerful combination of PFG NMR with magic angle spinning has recently 



enabled a notable increase in both observation times and gradient pulse intensities c 



As a 



consequence PFG NMR diffusion measurements became possible beyond the limits of mea- 
surability existing so far. This concerns in particular the first diffusion measurements with 
liquid crystals confined to nanopores To stimulate a direct experimental verification of 
our present predictions using these novel techniques is the primary purpose of this study. 
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